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In the past five years enormous progress has been made in the ab initio calculations of the optical 
■response of electrons in semiconductors. The calculations include the Coulomb interaction 
>- .between the excited electron and the hole left behind, as well as local field effects. However, they 
^ [are performed under the assumption that the atoms occupy fixed equilibrium positions and do 
'not include effects of the interaction of the lattice vibrations with the electronic states (electron- 
phonon interaction). This interaction shifts and broadens the energies at which structure in 
2 the optical spectra is observed, the corresponding shifts being of the order of the accuracy 
O claimed for the ab initio calculations. These shifts and broadenings can be calculated with 
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various degrees of reliability using a number of semiempirical and ab initio techniques, but no 
full calculations of the optical spectra including electron-phonon interaction are available to 
'date. 

This article discusses experimental and theoretical aspects of the renormalization of op- 



O .tical response functions by electron-phonon interaction, including both temperature and iso- 
[topic mass effects. Some of the theoretical techniques used can also be applied to analyze the 
renormalization of other response functions, such as the phonon spectral functions, the lattice 
^ parameters, and the elastic constants. 

1 Introduction 

The temperature dependence of macroscopic properties of crystals (i.e., specific heat 
thermal expansion thermal conductivity elastic constants [0) played an important 

role in the early development of the quantum theory of solids More recently, the tem- 

perature dependence of microscopic properties, involving spectroscopic elementary excitations 
such as those responsible for optical absorption and refiection and phonons [jlO| have been 



discussed both experimentally and theoretically. The theoretical treatments have ranged from 
semiempirical ones to state-of-the-art ab initio calculations [|lO[|ll . 
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During the past five years considerable progress has been made in calculating ab initio the 
linear response functions, i.e., the complex dielectric functions of semiconductors and insulators 



12-141. These calculations, which include the Coulomb interaction between the excited electron 



and the hole left behind in the valence band (the so-called excitonic interaction) and local-field 



effects, represent an enormous improvement over the early semiempirical calculations [p!5| ,p!q]. 
They nevertheless implicitly assume that the constituent atoms occupy fixed positions and 
thus neglect the effects of the lattice vibrations on the electronic structure, i.e., the electron- 
phonon interaction. Such effects have been calculated for selected transition energies (at so- 
called interband critical points, CP) [p|,p|,p!T|, but calculations of the full spectrum of electronic 
interband transitions including electron phonon-interaction effects are yet to be performed. 

Considerable effort has been devoted to experimental studies of the effect of external pertur- 
bations on the dielectric function of solids. Among these perturbations we mention static strain 
TTIjrSII , electric fields ||T^, magnetic fields [^, and sublattice displacements corresponding to 



optical phonons with q = [^. The latter determine the frequency dependence of Raman 
spectra. The contribution of a crystal surface to the linear optical response has also been pro- 
fusely studied. Several articles in these proceedings deal with this subject, both experimentally 
and theoretically, the theoretical work being not as advanced as that mentioned above for the 
bulk. In any case it does not include ab initio the effects of the electron-phonon interaction 
either. Nor do the rather simplified calculations which have been performed for the effect of 
strain |[T7|j2^j23[] or electric fields on the dielectric function. 



The effect of "temperature" on optical spectra of crystals results from the interaction of 
the electronic states with the vibrations of the lattice (electron-phonon interaction), an effect 
which does not vanish even for T = 0. According to quantum mechanics, the vibrational 
amplitude does not vanish at T = 0: the so-called zero-point vibrational amplitude remains 
and is responsible for energy renormalizations up to 100 meV, of the order of the accuracy 



claimed for state-of-the-art ab initio calculations of critical point energies in such spectra |12 



Zero-point vibrational amplitudes should vanish, in a gedankenexperiment, when the ionic 
masses of a crystal become infinitely large. In the past 12 years, many crystals with variable 
isotopic compositions have become available. It has therefore become possible to measure 
the dielectric functions e{uj) of crystals with different isotopic masses, to extrapolate them to 



infinite values of M and thus obtain experimentally the bare values of transition energies |25 



23. It may sound paradoxical that one can obtain bare, i.e., unrenormalized, parameters from 



experiment since only renormalized parameters, including all self-energy corrections such as 
those due to electron-phonon interaction, are observable. However, as pointed out by Allen [p9[] , 
the procedure to obtain unrenormalized parameters just described is only an approximate one. 
It requires assumptions which lead to a simple dependence of the zero-point renormalization on 
isotopic mass and temperature, in particular the truncation of the electron-phonon perturbation 
theory series after the second-order terms. 

In this paper we describe the three contributions which lead to the temperature depen- 
dence of optical transition energies (including zero-point effects) in crystals, up to second-order 
perturbation theory. They are: 

i) The effect of the first-order electron-phonon interaction Hamiltonian to second-order in 
perturbation theory. These are the so-called Fan terms P). 

a) The effect of second-order electron-phonon interaction in first-order perturbation theory 
||30|1. These are Debye- Waller terms, sometimes also referred to as Yu-Brooks terms 0]. 
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in) The effect of thermal expansion, including zero-point renormalization of lattice parame- 
ters, coupled with the dependence of the transition energies on such parameters 0]. 



We shall discuss simple algebraic expressions which have been used to fit the measured 
dependences of transition energies (also called gaps or critical points, CP's) on temperature, 
together with the related effects of isotopic mass on CP's. The effect mentioned above under [i) 
is described in terms of complex self-energies E = + ^ Sj, which encompass not only energy 
shifts (Er) but also broadenings of the optical features (Lorentzian-HWHM = — Ej). The effects 
(ii) and [iii) correspond to real energy shifts of CP's or energy gaps. Emphasis is placed in 
this article on a method of obtaining zero-point renormalization through extrapolation to zero 
temperature of the linear part of the shift of CP energies with temperature usually found at high 
temperatures. Extension of these phonon renormalization to other more macroscopic physical 
properties, such as lattice parameters (needed also for the evaluation of the {iii) terms), elastic 
constants, and phonon frequencies, is also discussed. 

2 Renormalization of Electronic States through Electro n- 
Phonon Interaction 

The renormalization of excitation energies by lattice vibrations (or other perturbing agents 
such as disorder) is represented, in the many-body theory of solids, by a frequency (= energy 
for h = 1, a value of h which will be assumed implicitly throughout this paper) and wavevector 
dependent self-energy E(co',k)= Er(a;,k) + i Ej(a;,k). While the bare excitation energy ci;(k) is 
real, i.e., the excited "particles" with a well-defined real wavevector k have an infinite lifetime, 
the imaginary part of the self-energy adds an imaginary component to the excitation energy 
and the excited "particle" acquires a finite lifetime r for a real value of k, i.e., it becomes 
a "quasiparticle" . This lifetime of the corresponding electron-hole excitation is given by the 
uncert ainty-principle-like expression : 

r = (-2 E,)-^ , (1) 

the minus sign resulting from the convention chosen for the sign of Ej (Ej < 0). We are 
concerned here with frequency renormalizations induced by phonons of frequency f2(q) and 
harmonic vibrational amplitude u{fl,q). Because the time average of u is zero (in the harmonic 
approximation, perturbation terms of first-order in u{Q,q) vanish in E(ci;,k)). Such terms 
are responsible for a broad, incoherent (i.e., that does not correspond to a single value of k) 
background which will not concern us here, except in so far as it has been incorrectly interpreted 
sometimes to generate a finite value of Ej(co',k) pl| . 



2.1 Thermal expansion effects 

As already mentioned, the calculation of E(co',k) is broken up into two parts. One of them 
involves a perturbation calculation keeping the solid at constant volume (i.e., keeping the 
shape of the primitive cell and the Brillouin zone, BZ, constant). However, the experiments 
are not performed at constant volume, but at constant pressure. Hence, we must add to the 
perturbation results just mentioned the effect of thermal expansion at constant pressure (usually 
atmospheric, for practical purposes zero): 
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_ ( AV{T) \ _ ( AV{T) \ (a<,A 

where AV{T) = V{T) — Vq is the change in volume induced by the temperature T and a is 
the change in transition energy ujq induced by the change in volume, the so-called hydrostatic 
deformation potential of the CP under consideration. Usually one measures the dependence 
of ujq on pressure, in which case Auq must be obtained from the r.h.s. of (1) using the bulk 
modulus Bq and the pressure dependence of ujq represented by (duo/dp). 

The thermal expansion (AV"(T)/Vo)p results from anharmonic terms in the expansion of the 
crystal energy vs. atomic displacement, whereby only odd terms in this expansion contribute 
(or combinations of odd and even terms, with an odd number of odd terms): The even terms 
average to zero and do not change the equilibrium position of the atoms. We shall keep for our 
discussion only the lowest (third) order anharmonic terms which lead to a linear dependence of 
AV{T) on T at high temperature, i.e., to a T-independent thermal expansion coefficient in this 
T range. These terms are conveniently represented by the mode Griineisen parameter 7n(q) 
defined as: 

where q is the wavevector of a generic phonon and n its branch index. Using themodynamic 
relations, and minimizing the free energy at "zero pressure" and temperature T, one finds 



32.33 



V{T) = Vo + -^T ^"(^) 7n(q)(2nB + 1) (4) 
where is the Bose-Einstein statistical factor: 



ns{n,T) = ;e=^, (5) 

with T in the same units as Q. Note that Eq. (4) allows us to calculate not only the temperature 
dependence of the crystal volume, but also the zero temperature renormalization, by setting 
riB = p2|-p^. In elemental crystals this zero-point renormalization is proportional to M~2, 



where M is the average isotopic mass. 

In order to evaluate the thermal expansion contribution to Aujo (including zero-point effects) 
with Eq. (2) one can use either experimental values of {AV(T) /Vo) or those calculated 

with Eq. (4). The former are obtained rather accurately by means of X-ray diffraction. For 
the latter procedure one can use either experimental values of 7n(q) |^ or values calculated 
using ab initio band structure techniques p^j33[] . One should mention, at this point, that the 
7n(q), which according to (4) should have an average value of about +1, is negative (~ —1) for 
the TA phonons at the edge of the BZ of most germanium and zincblende-type semiconductors 
(exceptions: diamond and, curiously enough, Cul [p8| , p9l ). When these negative values are 



entered into Eq. (2) they lead to anomalous negative expansion coefficients at temperatures close 
to the frequencies of the corresponding phonons. At much lower temperatures, the expansion 
coefficients become usually (but not always [Q) positive again. This results in a very flat 
dependence of V{T) vs. T for T below the TA-phonon frequencies [0,^. The oscillations 
due to the sign reversals of 7ra(q) are easy to see in the expansion coefficient vs. T which 
corresponds to the derivative of V{T) vs. T. They are hardly noticeable in AV{T) except that 
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the behavior of AV{T) vs. T is very flat at low temperatures. It is therefore very difficult to 
see the corresponding effects in the contribution to Aujq given in Eq. (2), a fact which, however, 
makes easier the analysis of data on the temperature dependence of electronic CP energies. 



2.2 Electron-phonon interaction 

The direct electron-phonon interaction terms, defined for a temperature independent volume, 
are responsible for the "explicit effect" of temperature on the electron energies. The term 
"implicit effects" is applied to the thermal expansion contribution discussed in 2.1. The two 
terms displayed as Feynman diagrams in Fig. 1 contribute to the phonon renormalization of 
one-electron states to the same order in perturbation theory as that used in 2.1. Electronic 
excitations involve an excited electron in the conduction band and a hole left behind in the 
valence band. It is easy to show that the renormalization can be obtained as the difference 
of two sets of diagrams of the type in Fig. 1, the electron diagram minus the hole diagram. 



provided one neglects vertex terms which connect electron and hole lines ||4T[|. Diagram (i), 
representing the first-order perturbation Hamiltonian taken to second order in perturbation 
theory, leads to a complex self-energy whose imaginary part Sj represents the Lorentzian width 
(HWHM) of the renormalized state. This Ej only vanishes at the uppermost valence state 
and at the lowest conduction state, i.e., for the lowest absorption edge (which can be either 
direct, like in GaAs, CdS, CuCl and GaN, or indirect, like in Ge, Si, diamond and GaP). 
The absorption spectra corresponding to this edge can, therefore, be very sharp at low T [ ^ . 
Diagram {ii) represents the second-order electron-phonon perturbation Hamiltonian, i.e., the 
second derivatives of the Hamiltonian with respect to the vibrational displacements, taken 
only to first order in perturbation theory and thus leading to simple frequency shifts without 
imaginary energy components. It is not uncommon to find in the older literature frequency 
renormalizations corresponding to only one of the diagrams in Fig. 1. Since 1983, however, 
both components of the electron-phonon interaction have been usually evaluated and added 

iii- 




(i) (ii) 

Figure 1: Electron-phonon interaction terms that contribute to the renormalization of electronic 
states, (i) Fan terms corresponding to a complex self-energy, (ii) Debye- Waller terms leading 
to a real energy correction. 

For elemental crystals the electronic energy shift induced by the term [ii) can be written as 



where A£)]y{n,q) is a real coupling constant specific to the electronic state under consideration 
and M the mass of the element at hand. In a k-dependent pseudopotential representation. 
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Eq. (6) can be evaluated by multiplying the pseudopotential form factors by Debye- Waller fac- 
tors [0. This procedure reduces simply to performing band structure calculations with two sets 
of pseudopotentials, the bare one and that screened by temperature-dependent Debye- Waller 
factors. This is the reason why calculations are often performed, using Occam's razor, only for 
the (u) terms, neglecting (i) |4^, a procedure known nowadays to be unjustifiable. Actually, 



when using a pseudopotential or any kind of truncated Hamiltonian for the calculation, the 
separate real parts of {€) and (m) do not have an independent meaning as they depend on the 
type of truncation at hand. Only their sum is physically meaningful. 

The self-energy terms (i) of Fig. 1 can also be represented by an equation of the type of 
Eq. (6) with a complex coupling constant AsE{n,q) which corresponds to the typical second- 
order perturbation expressions, with squared matrix elements in the numerator and an energy 
denominator with an infinitesimal imaginary convergence term in it . Its numerical evaluation 
is considerably more complicated than that of the Debye- Waller terms (iz), since it involves 
not only a sum over all phonon states, but also one over intermediate electronic states. The 
early analytical calculation by Fan |p, while pioneering, is therefore highly unreliable.^ Fan 
suggested a procedure to relate the (i) to the {iii) terms by assuming parabolic band edges in 
order to perform analytically the integration of Eq. (6). Two parameters, corresponding to the 
deformation potentials of the conduction and valence band edges, are then introduced. Their 
difference determines the {Hi) terms whereas the properly weighted difference of their squares 
determines the {€) terms. The (u) terms are ignored. This procedure, if correct, allows the 
separate determination of the hydrostatic deformation potentials of the valence and conduction 
bands as done in Ref. [43]. However, as already mentioned, the results of this procedure are 
not trustworthy. 

Note the similarity between Eqs. (6) and (4). At temperatures higher than the average 
phonon frequency (sometimes called the Debye frequency), they become linear in T. This 
fact, which follows from our assumption of terms only up to second order in the corresponding 



perturbation expansion P,|9|J4^ is often corroborated by experiments up to temperatures some- 
what below the melting point. For T — 0, Eqs. (4) and (6) tend to finite contributions to the 
zero-point renormalization. The high-temperature expansion of Eq. (6) can be written as 



1 

M 



(7) 



where the brackets represent average values. A similar expression applies to the high-temperature 
limit of Eq. (4). Considering the fact that in elemental crystals all frequencies Q are propor- 
tional to M~2, we conclude from Eq. (7) that in the high-T limit Au is independent of M 
(classical limit). For T — > (quantum limit) 



Au oc M"^ . (8) 

Equations (7,8) apply to all three components (z, ii, iii) of Au. If we are able to observe 
a linear region in Auj{T) at large T, the extrapolation of this asymptotic linear behavior to 
T = enables us to estimate the zero-point renormalization of the corresponding critical point 
(CP) energy Alu (T = 0). 

Figure 2 displays the temperature dependence of the energy of the lowest (indirect) gap of 
germanium and the linear high-T asymptote extrapolated to T = p6| , ^ . This extrapolation 
leads to the estimate Auj (T = 0) = -53 meV for the zero-temperature phonon renormalization 

"Fan obtained reasonable agreement with experimental data available to him, in spite of not having considered 
the (ii) terms. This is not unusual in the case of theories constructed so as to explain extant experimental results. 
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of the gap energy at K. Using Eq. (8) we can estimate a difference of 2.2 meV between the 
corresponding gaps of ^^Ge and ''°Ge, which compares well with the measured value of 2.2 
meV [0. Figure 2 displays not only the experimental points for Eg{T) but also the calculated 
sum of the (i) and {ii) contributions (dashed line) and this sum plus the thermal expansion 
contribution {iii) (solid line) [|^. 
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Figure 2: Temperature dependence of the energy of the lowest indirect gap of germanium. The 
dashed line represents the effect of electron-phonon interaction. The thick line represents this 
effect plus the thermal expansion effect. The thin line represents the asymptotic extrapolation 
to the bare gap at T = 0. The points are experimental EH . 



We should recall that the direct gap Eq of Ge is not its lowest gap, there being the indirect 
Eg gap (r — L) just discussed about 150 meV below. Hence, the zero-point self-energy of the 
direct gap Eq should have an imaginary part which should be small because this gap is not far 
from being the lowest. From Fig. 1 in Ref. 48 we estimate — Sj ~ 0.5 meV at T ~ K, much 
smaller than = 70 meV |^ . 

We have, thus far, lumped all contributions (i), (ii), {iii) to S into a single one. In order 
to illustrate that the three contributions have the same asymptotic behaviors, we display in 
Fig. 3 the measured change in lattice parameter a^iT) of diamond with temperature (related 
to ^V{T) in a straightforward way) |^|. The asymptotic linear extrapolation leads to a zero- 
point renormalization Aao/cto = 0.374 which, using Eq. (8), allows us to predict a difference 
between the lattice constants of (larger) and ^^C equal to 1.5 x 10~^, in excellent agreement 
with the experimental results and ah initio calculations |]50| , |32| . Notice that the one-oscillator 
fit of Fig. 3 leads to an average phonon temperature Vt = 1460K, which compares well with the 



corresponding value for the specific heat (1450K 0). The three oscillator fits of ||3^ leads to a 
dominant frequency Q = 2137K, rather close to that of a two-oscillator fit of the specific heat 
(1884 K @,i). 

Concerning the dependence of physical parameters on isotopic mass, in this section we have 
confined ourselves to elemental crystals: only the mass of one element appears in Eqs. (7,8). 
In the case of binary or more complex crystals, that dependence is different for each one of 
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Figure 3: Lattice parameter of natural diamond vs. temperature. The points represent 
schematically the experimental data (up to 2000K). The thin line is the high temperature 
asymptote which extrapolates to the bare ao for T ^ (horizontal arrow). The thick curve 
represents a three-oscillator fit whereas the dashed curve (only plotted for T > 2000K) corre- 



sponds to a one-oscillator fit. See [37| and text. 



the constituent masses and Eq. (7) becomes considerably more complicated. The experimental 



5T[| and theoretical information concerning dependences on isotopic masses for compound 
semiconductors is rather scarce. It can, nevertheless, be very interesting. In the case of the 
copper halide CuCl, for instance, the lowest gap (i?o) increases with increasing Cl-mass (a trend 
which agrees with Fig. 2) whereas it anomalously decreases with increasing Cu-mass pT| , |52| . 



Approximate Expressions for the Dependence of CP 
Energies on Temperature 



A number of empirical |^3[ and semiempirical pl|-[56[1 analytical expressions have been proposed 
to represent the measured Aci;o(T) and the corresponding imaginary part of the self-energy. The 
most popular fully empirical expression, proposed by Varshni, [KH] is: 



Aa;o(T) 



(9) 



At high T, Eq. (9) becomes linear in T, in agreement with the discussion of Sect. 2. For T — 0, 
however, it becomes quadratic in T, leading to a stronger variation with T than observed 
experimentally (see Fig. 2 in Ref. 56 for GaAs). 

The semiempirical expressions usually capitalize on the statistical factors (2^^ + 1) which 
appear in Eqs. (4) and (6), replacing them by average values. In Ref. 54 one average value, 
corresponding to a single Einstein oscillator as used in the theory of the specific heat 
is employed. In Refs. 43 and 55 two oscillators are used, paralleling the improved theory 
of the specific heat of Nernst and Lindemann [^. Like in the case of the specific heats 
electronic self-energies due to electron-phonon interaction at low T deviate from the single 
oscillator behavior. Two oscillators [@,^ improve matters a bit but, ub implying an exponential 
dependence on T, cannot reproduce the behavior found experimentally for the specific heat 
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of insulators. It is easy to prove that Eqs. (4), as well as (6), also lead to a low temperature 
behavior proportional to T^. This is accomplished in Eq. (4) by replacing the sum by an integral 
over energy, which requires the introduction of the density of phonon states, proportional to 
at low T. Switching to the dimensionless variable of intergration ^ = Vt/T leads to the factor 
that we are looking for.|] In order to obtain the same factor with Eq. (6), we need some 
additional reasoning. We must consider the fact that translational invariance requires that the 
electro n-phonon coupling constant vanishes for — > 0: long wavelength acoustic phonons are 
equivalent to a uniform translation which does not couple to the electronic states. Hence, a 
factor of VP' appears in A^y/ and Ase- The same factor of is then obtained in Eq. (4) for 
T ^ 0. We thus see that Eq. (9), leading to Aco'o(T) oc for T ^ 0, cannot be correct 
either. Several attempts to remedy this problem, using empirical expressions, appear in the 
works of Passler ||56|. We illustrate them in Fig. 4 for the indirect gap of Si. The experimental 



points fall above the single oscillator fit of Ref. 54, as expected. Also as expected, the fit with 
Eq. (9) falls well above the single oscillator fit of Ref. 54. An excellent fit is obtained with the 
four-parameter expression: 
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Figure 4: Temperature dependence of the indirect gap of silicon Eg{T). The points are experi- 
mental, the solid line represents a fit with Eq. 10. Other details are given in the figure and in 
the text. Adapted from 



(10) 



for p = 2.33. Equation (10) leads to a linear T-dependence for T ^ 6, as it should. Obtaining 
a good fit with an empirical or semiempirical expression is very helpful for the precise deter- 
mination of the linear asymptotic behavior which, as discussed in Sect. 2, can be used to find 
the zero-point gap renormalization. 



''Note that ifjn(<l) undergoes the changes in sign mentioned above, the range of behavior for the thermal 
expansion at low temperatures will be rather small. However, as already mentioned, these changes in sign make 
the behavior of V(T) vs. T very flat at low T, thus improving the overall quality of single oscillator fits, although 
they cannot describe the resulting non-monotonic behavior of the expansion coefficient for T — > 0. 
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An excellent fit to Eg(T), such as that provided by Eq. (10) (with p = 2.33, 9 = 405.6 
K, a = 0.3176 meV ■ K~^), can be used to circumvent the problem of determining the 
high temperature asymptote, and thus the zero-point renormalization mentioned above. Using 
Eq. (10) we obtain for the T = renormalization: 

AEo{T = 0) = - ^ . (11) 

Replacing into Eq. (11) the value a = 0.3176 meV ■ K~^, and G = 405.6 K given in [48], we 
find AEo{T = 0) = —64 meV for the indirect gap of Si. From this value of the zero-point 
renormalization one can estimate with Eq. (8) the isotope shift of the indirect gap of silicon, 
between natural silicon (M = 28.0856) and isotopically pure ^^Si, to be 0.12 meV, in good 
agreement with the experimental results (0.094 meV) [^]. 

We show in Table I the values of the zero-point renormalization of the lowest gaps (plus 
the lowest direct gap of Ge) of several diamond and zincblende-type semiconductors, together 
with the corresponding values of dE/dM. The first column under AE{T = 0) and dE/dM 
represents experimental data obtained by the linear extrapolation technique and by isotopic 
substitution, respectively (an exception is the case of diamond, for which the calculated value of 
AE{T = 0) is given). The available experimental data for diamond do not allow a reliable 



extrapolation to T = 0, the experimental value of dEg/dM listed in Table I is from 
The numbers in the second column under AE{T = 0) (and dE/dM) were estimated from the 
"experimental" values of dE/dM (and AE{T = 0)) using the M~2 dependence, as discussed 
above. The values of AE{T = 0) are rather similar for most materials in the table, with the 
exception of those for diamond, which are very large, and the rather small ones listed for the 
copper hahdes (to be discussed in Sect. 4.2). We note that several electron-phonon coupling 
constants of diamond are known to be considerably larger than those for other zincblende-type 



semiconductors a fact which we conjecture it arises from the absence of p-like electrons in 



the core of the carbon atoms. Correspondingly, the measured dEg/dM is also rather large for 
diamond 1571. 



4 Lowest Direct Gap Eq of compound Semiconductors: 
Temperature and Isotope Effects 

Information on the temperature dependence of the lowest gaps of many semiconductors can 
be found in standard compilations of semiconductor properties (e.g., the Landoldt-Bornstein 
Tables and Ref. [47]). We shall discuss in this section a few cases for which isotopic data are 
also available. 



4.1 GaAs and ZnSe 



The temperature dependence of Eq for GaAs is rather similar to that for germanium, as expected 
from the similarity in their band structures and lattice dynamics. The parameters for the best 
fit with Eq. (10) are 6 = 226 K, a = 0.473 meV -K'^ p = 0.251 [|9[. Using these parameters 
one predicts with Eq. (11) AEo{T = 0) = —53 meV, a value very close to that given in Table 
I for Ge. Unfortunately, arsenic has only one stable isotope (^^As) so that only the effect 
related to the two isotopes of Ga, ^^Ga and ^^Ga, can be experimentally investigated. Common 
sense, and also pseudopotential calculations |5^, suggest that the effect of varying the Ga mass 
should be nearly the same as that of varying the As mass (per unit atomic mass). The effect 
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of the ^^Ga-''^Ga substitution, realized experimentally, can thus be predicted by dividing the 
prediction of Eq. (8) for germanium by 2. We thus predict for the derivative of Eq with respect 
to the Ga mass: 

OE 1 

° X 53 = 0.19 meV/amu (12) 



dMca 4 X 70 

This value is a factor of two smaller than the measured one (0.39 meV/amu) and also than 
the value obtained from the pseudopotential calculations (0.43 meV/amu) [^. The reason for 
this discrepancy is not known but it may lie on having taken too small a temperature range for 
the fit in Ref. [^]. Indeed, if we use the data for Eq{T) given in Fig. 2 of Ref. for GaAs 



we find AEq{T = 0) =-90 meV, a value which leads to dE / dMca, in rather good agreement 
with the experimental result (see Table I). 

For ZnSe, the linear extrapolation of Eq{T) gives a zero-point renormalization AEq{T = 
0) = —65 meV [Q. Pseudopotential calculations yield a stronger electron-phonon con- 
tribution to the isotope effect for Se {OEo/dMse = 0.26) than for Zn {dEo/dMzn = 0.19),[| 
although the difference is compensated, at least partially, by the corresponding thermal expan- 
sion effects The measured values of dE^/dUzn = 0.214 and dE^/dMse = 0.216 are 
indeed nearly equal (see |l5Tl , |60| and Table I). The value estimated with Eq. (8) assuming that 



both derivatives are equal is 0.24 meV, in reasonable agreement with the experimental results. 

Notice that the thermal expansion contribution AEth{T = 0) shown in Table I for GaAs 
is about 40% of the total gap renormalization AEq{T = 0) and has the same sign. In the 
case of GaP, however, the thermal expansion contribution to the indirect gap is only 4% of 
the total effect and has the opposite sign. This apparent anomaly is related to the large value 
and the negative sign of the deformation potential of direct Eg gaps and the small value plus 
sign reversal which takes place for the F — > A indirect gaps Eg. Similar effects can be seen in 
Table I for AEq{T = 0) of Ge as compared with AEg{T = 0) of Si. For the indirect gap of 
germanium Eg, the sign of the thermal expansion renormalization is the same as that of Eq 
and its absolute value is somewhat smaller (25%) as corresponds to a F — >• L indirect gap. 

4.2 The cuprous halides 

CuCl, CuBr and Cul (I- VII compounds) crystallize, under normal conditions, in the zincblende 
structure. However, when compared with their group IV, III-V and II- VI counterparts, they 
display a number of anomalous properties which are usually attributed to a strong 

hybridization of the halogen p valence bands with the 3d electrons of copper. Among these 
properties we mention, as relevant to the present paper, the increase in Eq with increasing T 
observed for CuCl and CuBr.0 



Figure 5 displays the anomalous temperature dependence of the Eq gap of CuCl [^. The 
experimental points show an increase in gap energy with increasing temperature, with an in- 
creasing slope at low T that begins to decrease at ~100 K and seems to saturate at ~300 K. This 
behavior can be phenomenologically interpreted by means of a two-oscillator fit 0, provided one 
assumes that the two oscillators give electron-phonon contributions of opposite signs: The lower 
frequency oscillator, most likely related to Cu vibrations, contributing a positive AEq{T = 0) 
and the high frequency one, related to CI, giving a negative contribution. These opposite signs 
agree with those determined experimentally for the derivatives of Eq with respect to the copper 



'^Note that this trend continues for CuBr where dEo/dMcu even becomes negative J5 
'^For Cul Eq decreases only very slightly between K and room temperature poj^' 



11 




100 200 300 
Temperature [K] 



Figure 5: Temperature depen- 
dence of the Z3 exciton energy 
(equivalent to that of the Eq di- 
rect gap) of CuCl. The points 
are experimental, the thick solid 
curve represents a two-oscillator 
fit (acoustic and optic) while 
the thin line indicates the linear 
asymptote. From . 



and chlorine masses [dEo/dMcu 
Fig. 6.) PI. 



-0.076 meV/amu, OE^/dMci = -1-0.36 meV/amu, see 
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Figure 6: Two-photon absorption 
in the region of the longitudi- 
nal exciton (LEI) in isotopically 
modified CuCl measured at 2K. 
The isotopic shifts shown in this 
figure are expected to correspond 
to those of the direct gap Eq+Aq, 
which should be nearly the same 
as those of Eq. Notice that the 
shifts with increasing mass are 
negative (i.e., anomalous) for Cu 
and positive for CI 



The question arises as to whether the anomaly just mentioned is due to the electron- 
phonon interaction or to the thermal expansion term. Fortunately, the deformation potential 
dEo/dlnV , and also the thermal expansion, are very small and the corresponding effect can be 
neglected for CuCl and approximately neglected for CuBr, see Table I. The anomalous contri- 
bution of Cu, and that of the halogen, must then be due mainly to electron-phonon interaction. 

The two-oscillator contributions to AEo{T) are labeled "acoustic" and "optic" in Fig. 5. 
The anomalous acoustic contribution has been attributed in Ref. 51 predominantly to the Cu 
vibrations, the optic contribution, with the usual sign, to chlorine. An excellent fit to the ex- 
perimental data is obtained when adding both contributions. The horizontal line represents the 
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bare gap as obtained from the two-oscillator fit. A thin straight line represents the asymptotic 
behavior at high T and extrapolates to the bare gap for T = 0. 

A similar anomalous behavior of Eq{T) and OEq/Mcu is found for CuBr pll , |65| , |66| . A 
negative dEo/dMcu is also found for CuI.Q However, at least below 300 K for Cul, Eq{T) 
decreases with increasing T . This case is in the process of being analyzed . 



5 Higher Critical Points: Dependence on Temperature 
and Isotopic Mass 

We have labeled as Eq the lowest direct critical point (i.e., energy gap) of a germanium and 
zincblende-type semiconductor. In all cases discussed, Eq is found around the F-point (center) 
of the Brillouin zone. If it is the absolute minimum of all transition energies, including indirect 
ones (e.g., for GaAs, see Fig. 7, corresponding Ej vanishes.^ If not, a finite value of 

Ej is expected (e.g., Ge, — S,; = 0.5 meV at T = 0, see Sect. 3). The Eq CP at F is split by 
spin-orbit interaction, the second component being known as Eq + Aq. A finite value of Sj is 
expected for this second component since it is not the lowest: it is induced by the electron- 
phonon interaction between the coresponding valence state at the F point and the top valence 
band of F, at nearly the same energy p7|-[70| {{i) terms of Fig. 1). In the case of GaAs, — Ej 
for Eq + Aq has been found to be 6 meV at 2K [^,^. The behavior of -Sj with increasing 
temperature can be estimated by multiplying 6 meV by 2(72^) + 1, (jib) being the Bose-Einstein 
factor for the appropriate average phonon frequency. 
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Figure 7: Spectra of the dielectric function e2{uj) of GaAs ellipso metrically measured at several 



temperatures. Notice the critical points labeled i^i + Ai, Eq, and E2. From ||6^ 



Eq and Eq+/S.q correspond to optical absorption coefficients of the order of 10^ cm~^ (e2 — 1). 
With increasing photon energies the absorption coefficient rises to a maximum near 10^ cm~^: 
The corresponding CP is labeled E2. The €2(0;) spectrum of GaAs in the region of strong 
electronic absorption (1.5 to 5.5 eV) is shown in Fig. 7 for several temperatures. The strongest 
feature corresponds to the E2 CP's: it red-shifts and broadens strongly with increasing T. 
Another feature, labeled Eq, corresponds to transitions between the F25' valence band and the 
Fi5 conduction band states [|67|. 



e of natural (i.e., isotopically mixed) Si has 
d thus can be attributed to isotopic disorder 



'^Unfortunately, iodine has only one stable isotope (^^^I). 
fNote, however, that a broadening of the indirect absorption edg 
been recently discovered. It disappears for isotopically pure "^^ Si an 
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We display in Fig. 8 €2(0;) spectra measured for GaAs at 20 and 300 K (from Fig. 7) and 
compare them with recent ab initio calculations of Rohlfing and Louie ||12| , one of the calculated 
spectra (dashed line) without electron-hole (exciton) interaction. The measured spectra display 
clearly the spin-orbit splitting of the Ei transitions {Ei — i?i + Ai). The spin-orbit interaction 
was not included in the calculations in order to reduce (by a factor of 2 x 2 = 4) the size 
of the Hamiltonian matrix. The calculations include, however, an artificial broadening F ^ 
0.10 eV so as to smooth out the limited k-space integration mesh. This is the reason why the 
calculated Ei — Ei+ Ai structure, in principle not including lattice vibrations, is broader than 
that measured at 20K. Unfortunately, because of the single isotope nature of arsenic (^^As) and 
the rather small mass difference between the two isotopes of Ga (^^Ga, ^^Ga), the experimental 
accuracy does not allow for an extrapolation of the e2{uj) data to infinite isotopic masses, a 
procedure which would be desirable so as to compare measured with calculated spectra. In the 
case of Ge, however, the isotopic mass range available is much larger (™Ge— ™Ge). Moreover, 
one gains an additional factor of 2 (with respect to GaAs) because for Ge both atoms in the 
primitive cell can be isotopically substituted. Nevertheless, such an isotopic extrapolation 
procedure, which we believe to be possible but difficult, has yet to be performed. 




0) (eV) 



Figure 8: Spectra from Fig. 7 at 20K (blue) and 300K (red) compared with the calculations of 
Rohlfing and Louie [0 with (thick solid line) and without (dashed line) excitonic interactions 
but without spin-orbit coupling. The calculations include an ad hoc Lorentzian broadening 
of 0.15 eV. Notice the enhancement of the Ei and Ei + Ai peaks induced by the excitonic 
interaction. 



We shall discuss here in some detail the temperature and isotope effects of the Ei and 
El + Ai features of Ge and Si. They correspond to transitions between the top valence bands 
and the bottom conduction band along the (111) directions (labeled as A). Their splitting, 
190 meV for Ge, is due to spin-orbit interaction. Similar features are observed for all materials 
of the germanium-zincblende family |^|: A red shift (Er) and a broadening with increasing 
temperature. 
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5.1 The El and Ei + Ai critical points of germanium: temperature 
and isotope effects 

Ronnow et al. were able to determine with a remarkable accuracy the dependence of 

S(T = 0) on isotopic mass for the Ei and Ei + Ai CP's of germanium. We show in Fig. 9 the 
third derivative of ellipsometrically determined spectra of ei{uj) and £2(0;) for three germanium 
wafers with different isotopic compositions measured at 30K. It is easy to see in this figure that 
the El and E'l + Ai curves shift to larger frequencies with increasing isotopic mass, in agreement 
with the results found for the Eg and £^0 critical points (Table I). The critical point parameters 
El, El + Al, and F (HWHM) obtained from the spectra of Fig. 9 are plotted in Figs. 10 and 
11 versus isotopic mass. In spite of the considerable scatter of the experimental points, a clear 
trend is apparent: Ei,Ei + Ai, and even their difference Ai, increase with increasing mass. 
The solid lines in Fig. 10 are fits to E = E^o + BM^ , i.e., to the prediction of Eq. (8). These 
fits enable us to determine the bare values of the critical points E^o] they are listed in Table II. 
Notice that the spin-orbit splitting Ai shows a slight tendency to increase with increasing mass, 
surely a result of the electron-phonon interaction; quantitative calculations of this increase are 
not available. 
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Figure 9: Third derivative with 
respect to photon energy of the 
dielectric function of natural Ge 
in the vicinity of the Ei and Ei + 
Al transitions. The solid lines 
are fitted to the experimental 
data using analytical line shapes. 
Red curves: ™Ge, black; natural 
Ge, green; ^^-^Ge, |71. 
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Figure 10: Energies of the Ei and Ei+Ai 
transitions in Ge, vs. isotope mass. Also 
shown is the dependence of the spin-orbit 
splitting, Al, on isotope mass. The solid 
lines are fits io E = E^+BM'^ (see W\\). 



Figure 11 displays the decrease in the Lorentzian HWHM parameters of the Ei and Ei + Ai 
critical points of germanium observed with increasing isotopic mass. The solid lines are fits 
with the equation F = Foo + BM~2, based also on Eq. (8), except that a small intercept Fqo 
(^!^^ —6 meV) allowed us to improve the fit. This intercept, if meaningful, would represent a 
"natural width" which would be present even when the atoms do not vibrate. It could be due 
either to crystal imperfections, among them the truncation at the surface, or, alternatively, 
to inadequacies in the line shape used to fit the experimental spectra. In Ref. 71 logarithmic 
singularities were used (see Eq. 4 in [^) which correspond to two-dimensional interband critical 
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points. Excitonic effects may contribute a small negative Too to the fitted r(M) as also may 
the lowered dimensionality induced by the surface. The fitted value of B is 350 ±90 meV/amu, 
in remarkably good agreement with band structure based calculations (310 meV/amu [@]). 




isotope mass [amu] 

Figure 11: Lorentzian widths F = — Sj of the Ei and Ei -\- A transitions of germanium vs. 
isotopic mass. The points are experimental. The lines represent fits to F = Too + BM~2 [[71|] . 



5.2 The El critical points of silicon: temperature and isotope effects 

Work similar to that described in Sect. 5.1 for germanium has been performed for silicon by 
Lastras, et al. using ^^Si, ^°Si and natural silicon (M = 28.09). Because of the small effects 



involved, the fits to the experimental lineshapes were performed with an elaborate, highly 
accurate Fourier transform technique |]72|. The analysis of the data is further complicated 
by the inability to resolve the spin-orbit splitting Ai (for 5'j,Ai 30 meV), and the near- 
degeneracy of the El CP's with the critical point Eq. Therefore, it was not possible to obtain 
data for the dependence of Sj on isotopic mass. The width of the Ei transitions measured 
directly at T = 10 K is about 70 meV |7^, a value somewhat larger than that found for the 



other materials under discussion. The value of Er(T = 0) obtained from the isotope shifts is 
118±22 meV, somewhat larger than the calculated value (71 meV) and much larger than the 
value extracted from the measured temperature dependence of the Ei transition (39 ± 13 meV) 
73|. The discrepancy may be due to the complex nature of these transitions, as mentioned 



above. 

We display in Fig. 12 the results of a recent ah initio calculation of e2{oj) for silicon, per- 
formed both with (solid black line) and without (dashed line) excitonic interaction [p!^ 



com- 



pared with spectra measured ellipsometrically at 20 and 300K |7^. This figure emphasizes 



the importance of comparing calculations performed for static atoms with measurements at 
the lowest possible temperatures: the height of the Ei peak measured at 20K agrees rather 
well with the calculated one. The calculated peak is, however, broader than the measured one. 
This is a trivial effect of having introduced in the calculation an artificial Lorentzian broadening 
parameter of about 100 meV. The calculated Ei structure is of obvious excitonic nature since 
it goes down dramatically in strength if excitonic interaction is not included. It peaks 100 
meV higher than that measured at 20K. While 100 meV is about the accuracy claimed for ah 
initio calculations, it is tempting to speculate, in view of the results for Ai?i(T = 0) given in 
Table II, that the red shift of 100 meV found in the experimental spectrum with respect to the 
calculated one is due to the zero-point renormalization by the electron-phonon interaction. 
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Figure 12: Ellipsometric spectra of e2{uj) for silicon [^] at two temperatures compared with the 
calculations of Albrecht et al. |]T3 with (solid black line) and without (dashed line) excitonic 



interaction. Note that the peak calculated at 3.8 eV does not appear in the experimental data 
and thus it is to be regarded as spurious, possibly due to the coarse a k-space interaction mesh. 



There has been some controversy concerning the peak exhibited by the black solid curve of 
Fig. 12 at 3.8 eV. While in the original publication |jl^ it was claimed that this peak corresponds 
to a well-established experimental feature, Fig. 12 contradicts this claim: It seems that the 3.8 



eV feature is simply an artifact of the coarse mesh used for the calculation [25,75 



6 Linear Optical Response Related to Surface States: 
Effects of Electron-Phonon Interaction 

There have been in recent years a large number of investigations of the response of electronic 
surface states to near-ir and visible light. This response is rather weak and can best be mea- 
sured in a differential way which involves either subtracting the responses to two inequivalently 
polarized light beams, or letting the surface states disappear through passivation (e.g., oxida- 
tion 1^). Even so, measurements in uhv are difficult and time consuming. Consequently, only 



few experimental investigations of the temperature dependence of the corresponding spectra 



have been published |3T|J77|J78[] . These investigations are mainly concerned with the gap be- 
tween filled and empty electronic surface states of cleaved [111] silicon surfaces. For a critical 
discussion see Ref. 79. 

The spectra of 2 x 1 [111] silicon surfaces cleaved in vacuum consist of a strongly polarized 
(for E[[ [110], no response for E[[ [211]) peak centered at about 0.5 eV. This peak barely shifts 
with temperature, a fact which has been attributed in Ref. 79 to a compensation of the shifts 
induced by the various phonons involved.^ 

A red shift is, however, found for the low-frequency onset of this peak (0.43 eV at 20 K, 0.34 
eV at 300 K |7^). The magnitude of this shift is similar to that found for the lowest absorption 
edge of many semiconductors (see Figs. 2 and 4). The onset of a weak absorption peak, however, 
cannot be very reliably measured and considerable differences exist among values reported for 



^This is similar to the situation found for bulk Cul 
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Si. Nevertheless, several attempts have been made to interpret the red shifts reported for the 
onset and the near lack of shift at the peak |^ . The authors of these attempts basically accept 
the latter fact and explain the shift of the onset as a broadening of the peak structure. 

An interesting observation concerning the response of the 2 x 1 [111] silicon surface is the 
broadening with increasing temperature which can, of course, be attributed to electron-phonon 
interaction. We have seen in the previous sections that this broadening is proportional to T at 
high temperatures and becomes T- independent for T — 0, at least for bulk semiconductors. In 
the case of surfaces one may argue that, because of the reduced dimensionality, the coupling to 
the surface electrons of individual surface phonons may be much stronger than in the bulk and 
therefore the lowest order perturbation theory which led to Eq. 6 may break down ||8^. For 
the [111] silicon surface there is no convincing evidence for this to be the case. 

It has been mentioned in Sect. 2.2 that, to lowest order in perturbation theory, both AE{T) 
and r(T) are proportional to (w^) {u = phonon amplitude), i.e., to 2{nB) + 1, where (ub) is an 
average Bose- Einstein factor (Eq. (5)). This fact leads immediately to the linear dependence 
on T at high T which has been so profusely and fruitfully used throughout this article, the 
surface spectroscopist will, however, have found in the literature statements to the fact that 
the widths of spectral structures induced by transitions between electronic surface states are 
proportional to [2(nB) + l]^/^, i.e., to T^/^ at high T [31,79[. This is a point where communication 
between surface and bulk physicists could have been better. The argument behind the linear-w 
dependence emphasizes the fact that the vibration, symmetric with respect to u = 0, slushes 
back and forth the frequency of the peak, thus producing a broadening proportional to u. Those 
familiar with the Debye- Waller theory will, however, realize that such a term generates only 
a broad background, not a renormalized band-structure preserving sharp critical points. The 
renormalization of the band structure contains, in the spirit of Fig. 1, only terms proportional 
to 2{nB) + 1. 
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Figure 13: Widths (HWHM) of the spectral peak observed around 0.5 eV for a 2 x 1-Si [111] 
surface at four temperatures. The three points above 100 K are from [^, that near T = was 
recently communicated to us by the authors of [pT|. The line through the points represents a 



fit with the function [2nB{Q/T) + 1] x r(0) which led to 6 
for silicon. 



340 K, a rather reasonable value 



Unfortunately, the only measurements vs. temperature of the width of the 2 x 1 Si [111] 
peak were reported solely for three temperatures [^ (see Fig. 13; a fourth point, close to T = 0, 
has been recently communicated to us by the authors of [^) and they hardly possess enough 
information to refute experimentally the square- root hypothesis. In fact, they were presented 
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in |31] as supporting it. We have plotted in Fig. 13 these points vs. T and fitted them with 
a single-oscillator 2nB + 1 expression. The fit for an average frequency 6 = 340K, is quite 
satisfactory, slightly better than the square-root fit of Fig. 2 of |3l[ . 



7 Dependence of Other Physical Properties on Temper- 
ature and Isotopic Mass 

The Bose-Einstein fits and measured temperature dependences of a wide variety of physical 
properties of semiconductors, in connection with linear asymptotic expansions, can also be 
used to estimate the zero-point renormalization of these properties and the corresponding iso- 
tope effects. In this section we briefly discuss the long-wavelength refractive index, the phonon 
frequencies and the elastic constants for elemental semiconductors. The semiempirical expres- 
sions used for the fits are basically the same as those discussed in previous sections for interband 
transition energies. However, the M~2 dependence of the zero-point renormalizations, given 
in Eq. (8), must be modified depending on the type of property under consideration. This 
dependence still holds for the long wavelength refractive index since it is directly related to 
ei{uj) and e2{uj) (a straightforward calculation suffices to prove it). 

In the case of phonon properties, such as phonon frequencies, the M~2 dependence must 
be replaced by M~^. This is easily understood. The diagram, corresponding to (ii), which 
represents an anharmonic renormalization of a phonon frequency, is obtained by replacing the 
horizontal electron propagators by curly phonon propagators. Hence, four phonon amplitudes 
come together at the vertex, giving a matrix element proportional to u^, i.e., to M~^. For the 
elastic constants, proportionality to M~2 obtains. 



7.1 Long wavelength refractive index: temperature dependence and 
isotope effects 

The temperature dependence of the long wavelength refractive index plays an important role 
in the tuning of optical and optoelectronic devices. Early theoretical work was concerned only 
with the temperature region in which the T-dependence is linear P5|j81| . Recently, ab initio 
calculations over a broad temperature range, including the zero-temperature renormalization, 
have been performed Experimental results are given in for diamond and in |^ for Si 
and Ge. 

The following fitted expression is given in Ref. 82 for the long wavelength refractive index 
n{T) of diamond. 



n(T) 



A 

no+ - 



2nB (fi,T) 



(13) 



where Q = 711 cm~^, A = 0.19 and uq = 2.377. This expression implies that the zero- 
point renormalization is An(T = 0) = 0.085, rather close to the value 0.07 obtained by linear 
extrapolation from Fig. 5 of Ref. 7. Note, however, that the value of An given in Table I of 
Ref. 7 is three times larger and thus must be erroneous, since inconsistent with Fig. 5 of that 
reference. 

From the zero-point renormalization An(T = 0) = -1-0.085, the dependence on isotopic 
mass: 
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dno 
dM 



4 X 10"^ K^' (14) 



is obtained. It is slightly larger than the value 3 x 10~^ calculated in Ref. 7. Fabry-Perot 
interference fringes are a very accurate method to measure differences in refractive indices. 
They could be used to measure the coefficient of Eq. (14) by comparing two plane-parallel 
diamond wafers made out of ^^C and ^^C, respectively. The expected value of this coefficient, 
however, is so small that differences in the samples resulting from their fabrication, but due to 
reasons other than their isotopic mass, would be likely to obliterate the expected results. At 
this time it would also be very difficult to measure the corresponding derivatives of uq for Ge 
and Si which, nevertheless, can be easily obtained from the calculations in [0. 

7.2 Anharmonic self-energy of phonons: effects of temperature and 
isotopic substitution 

Although anharmonic self-energies can be found for all phonons in the dispersion relations, 
most of the information available concerns the Raman phonons, i.e., the optical phonons for 
q = 0. In most cases analyzed in detail, the equivalent of the two terms in Fig. 1 (with the 
electron lines replaced by phonon propagators) suffice to explain the measured T,{Q) [p!0| , p^ ,^ . 

Let us first discuss the dependence of the phonon linewidth (HWHM = — Sj) of the equiv- 
alent of Fig. 1 for phonons. No contribution of the phonon-phonon equivalent of (ii) exists 
|]85|j86[| . The contribution of (z) to Sj represents a decay of the phonon under consideration 
(frequency fl, q=0 for the Raman phonon) into two phonons, with frequency and wavevector 
conservation {Q = Q1 + Q2] q = ^1 + ^2)- In the case of Ge and Si excellent fits to the measured 



temperature dependence of the HWHM are obtained with the expression p5|,|86[]: 



r(T) = -S,(T) = Fo [1 + nB{n^,T) + 71^(^2, T)J (15) 

where Qi = 2Q2 = 0.66^7 {Qi and ^2 correspond to LA and TA phonons at the edge of the 
BZ), and the zero-point renormalization Fq — 0.5 cm~^, which has indeed been shown to be 
proportional to M~^.Q The energy conservation constraint i7 = + ^2 is only required for 
the evaluation of E,. 

The frequency shift S^, including the contribution equivalent to (ii), is obtained by summing 
over all possible virtual intermediate states, with (q = qi -l- q2) but without the requirement 
of energy conservation. However, because of the energy denominators which appear in the 
second-order perturbation expression corresponding to S^., two phonons with Qi + ^2 near Q 
often give the dominant contributions to S^. Hence, it is common to find in the literature fits 
to the measured Sr(T) with Eq. (15) and the restriction = fix -|- Such a fit is found 
for silicon in Fig. 3 of Ref. 83. It leads to an estimate = —5.6 cm~^ for the zero-point 
anharmonic renormalization of the phonon frequency of Si (bare frequency at T = equal to 
531 cm~^). This zero-point renormalization, proportional to M~^, decreases in absolute value 
by 0.39 cm~^ when replacing ^'^Si by ^°Si. 



'^See Fig. 4 in Ref. 83, where the contribution of isotopic disorder to T has also been included. This contri- 
bution amounts to 0.035 cm~^ HWHM for a mixture of equal parts of'^^Si and '^'^ Si. 
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7.3 Elastic Constants 



The temperature dependence of the elastic stiffness moduh Cii(T), Ci2(T) and Cm{T) has 
been measured for diamond sihcon ||88|, and germanium |8^. In spite of the temperature 



range, which was hmited to T shghtly below the Debye temperature, it has been possible to 
estimate the position of the linear asymptote required for the determination of the zero-point 
renormalization. As already mentioned, this renormalization is related to that of the long- 
wavelength acoustic phonons and turns out to be proportional to M~2. 

In Table I of [49] the values of the elastic constants, the estimated temperature renor- 
malization and their logarithmic derivatives with respect to M, calculated by using the M~2 
dependence on mass, are displayed. No experimental data are available for these derivatives 
but, in spite of the small values involved, they may be measurable using the highly accurate 
ultrasonic propagation technique. 

It is interesting to mention that a value of dlnC/dM = 5 x 10~^ was reported in for 
diamond. A point was made in that paper of the fact that correspondingly, ^^C diamond must 
be significantly harder than ^^C diamond. This dlnC/dM is about an order of magnitude 
larger than the value obtained from Table I of Ref. [49], which implies insignificantly harder 
^^C diamond. In a following paper [90] an error in the interpretation of the experiments of 
|89[| was corrected. This correction, and an estimate with a theoretical model, led to a value of 
dlnC/dM < 10~^, in agreement with the value given in Table I (3.5 x 10"'^). 



8 Conclusions 

We have discussed the renormalization of a broad range of physical properties of tetrahedral 
semiconductors by the electron-phonon interaction. These properties include the dielectric func- 
tions and the corresponding electronic transition energies (both bulk and surface contributions), 
the lattice parameter, the phonon frequencies and the elastic stiffness constants. Particular at- 
tention has been paid to the zero-point (T = 0) renormalizations, which can be changed by 
changing the isotopic masses. Semiempirical expressions representing these renormalizations 
vs. temperature have been given. Especially powerful for the interpretation of experimental 
results has been shown to be the determination of the asymptotic linear behavior vs. tempera- 
ture at high temperature, either with the naked eye or, even better, from a fit to a well-founded 
analytical expression. 

In this article we have only discussed elemental and binary semiconductors. A few cases 
involving materials with three different elements or more can be found in the literature. Of 
particular interest are the I-III-VI chalcopyrites containing as element either copper or silver 



(e.g., CuGaSe2 which is isoelectronic to ZnSe; AgInSe2 ||9^ isoelectronic to CdSe; AgGaSe2 
|93[| ) . These materials all exhibit anomalies in the temperature dependences of the gap which are 
reminiscent of those observed for the copper halides: at low T the gaps increase with increasing 
T, whereas they decrease at higher T's. While it is easy to conjecture that the anomalies are 
due to the presence of 3d electrons of Cu (or 4d of Ag) in the valence bands, measurements 
of the dependence of the low temperature gap on isotopic mass would be desirable in order to 
tune these anomalies. Many of these chalcopyrites (e.g., Cu,Ga,Se) have three elements with 
more than one stable isotope each. 

The present author was delighted to see, after completing this paper, a very recent article 



94(1 in which the temperature renormalization of eight prominent structures, seen from 2 to 5 
eV, in the reflectance anisotropy spectra of (100) InP surfaces were reported. The experimental 
data for the temperature shifts were fitted with single oscillators, i.e., with expressions similar 
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to Eq.(6). The zero-point renormalizations of most of the eight observed structures he around 
-50 meV, rather similar to those of the bulk CP's. The fitted oscillator frequencies are close 
to the average phonon frequencies of bulk InP. As required by the oscillator fit, all observed 
frequencies depend linearly on T at high T. A square- root dependence is thus ruled out, at 
least for the frequency shifts. 
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Table I: Zero-point anharmonic renormalization AE{T = 0) of the lowest gaps of several 
zincblende-type semiconductors (in meV) together with its derivative with respect to the iso- 
topic mass (in meV/amu). Direct and indirect gaps are labeled Eq and Eg, respectively. For 
completeness, the thermal expansion contribution to the zero-point parameters is also listed 
{AEth and dEth/dM). For the binary compounds the derivatives with respect to cation (above) 
and anion (below) masses are given. The first number under AE{T = 0) has been measured 
through linear extrapolation, except for diamond where a value calculated in is given. 
The first number under dE/dM has been measured through isotopic substitution. The second 
numbers under AE and dE/dM were estimated from the measured values of dE/dM and AE, 
respectively, using Eq.(8). For sources see text. 
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Table 11: Derivatives of the Ei and Ei + Ai critical point energies of germanium at low tempera- 
tures as obtained from Figs. (9) and (10). Also, corresponding zero-point gap renormalizations 
estimated from the above values with Eq.(8) (in meV/amu and meV) and corresponding values 
for the El critical point of silicon |2^. The experimental values measured by isotopic substi- 
tution are nearly a factor of two larger than those calculated using empirical pseudopotentials 
MM. 
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